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DETERMINATION OF SATELLITE ORBITS 
FROM RADAR DATA 


SUMMARY 

An optimum method for determining satellite orbits from radar data is presented in 
this report. Offering a good combination of speed and accuracy, the method makes 
use of orbit inclination and orbit elements in the plane, and proceeds with a differential 
correction of the orbit elements. 

Rapid, accurate methods of computing orbit elements are required to predict satel- 
lite positions for acquisition by other radars at points later along the orbit. In some cases 
the data are limited to a single pass over the observing station. The dynamical method 
is described in detail, and its accuracy is compared with those of two other methods: 
the purely geometrical, and the least- squares geometrical. 

By this optimum method the computing time, including the differential correction 
time, is 1 minute. Without differential correction, the rough determination takes from 
10 to 20 seconds with approximately 5 miles positional uncertainty. 
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DETERMINATION OF SATELLITE ORBITS 
FROM RADAR DATA* 


INTRODUCTION 

An investigation of methods for determining satellite orbits from radar is described 
in this paper. Particular emphasis is placed on problems that require rapid calculation 
of the orbit elements from limited data. In some cases the available information may be 
restricted to the data obtained from a single pass, over a period of 2 to 4 minutes, as the 
satellite goes by the observing station. The solution to this computing problem requires 
rather different methods than those previously developed for calculating orbits based on 
several satellite revolutions. 

A procedure which offers a good combination of speed and accuracy is described in 
some detail in the following pages. This method requires between 10 and 20 seconds for 
a rough determination, and yields orbit elements in that time with an accuracy of 0.001 
in eccentricity and 0.05 minute in period— a degree of accuracy sufficient to determine 
whether the satellite is in orbit. The program then proceeds automatically with a 
correction routine that operates on the preliminary orbit elements to yield a more 
accurate prediction of the satellite position at later points along the orbit. The accuracy 
of these corrected predictions is 500 yards or better, which is sufficient for acquisition 
by other radars at later points along the orbit. This correction routine requires an 
additional 40 seconds of computer time. 

PROCEDURE 

The procedure chosen as optimum involves the following steps: 

(1) Orbit inclination 

A plane is passed through the center of the earth and through all the points along the 
radar track. The inclination of the plane and the position of the line of nodes are adjusted 
to the full set of data by a least-squares calculation. The plane of the orbit is then 
determined. 

(2) Orbit elements in the plane 

After the plane of the orbit has been fixed, the three orbit elements in the plane must 
be determined. These may be taken as the period, perigee altitude, and argument of 
perigee. They can be computed from any three items of information along the track. 

The most accurate way found was to choose two points on the track and the time between 

^Originally presented at the First International Space Symposium sponsored by the Com- 
mittee on Space Research (COSPAR), Nice, France, January I960. 
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An IBM-704 time of 1 minute is required for the determination of the orbit elements 
from one radar pass including twenty data points. By elimination of the differential cor- 
rection routine, the computing time can be reduced to between 10 and 20 seconds at the 
sacrifice of a factor of 10 in the accuracy of the orbit determination. In this case, the 
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them as the three fundamental items. The procedure for calculating the orbit ele- 
ments depends in an essential way on the dynamics of the satellite motion and may be 
called the "dynamical" method. 

f..n io rnorlo of fho rlata hv riiviHinc the track into halv es and nairintr the DOints 
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DETAILS OF THE ANALYSIS 

The x , y , z Cartesian coordinates of the satellite are determined from the slant 
range, bearing, and elevation. The orbit plane is then determined by a least-squares 
adjustment of the angle of inclination and the position of the line of nodes at a given time. 
Appendix A gives the equations used for the reduction of the radar data. 

The second step in the procedure consists in the determination of the planar orbit 
elements, that is, the period, eccentricity, and argument of perigee. The dynamical, or 
Gauss-Olbers, method has been found to give the most accurate results in this deter- 
mination. 


GAUSS-OLBERS METHOD FOR DETERMINATION OF THE PLANAR ELEMENTS 

The equation of an ellipse expressed in terms of the eccentric anomaly is 


r = a( 1 - e cos E) , 


and the Keplerian equation is 


( 1 ) 


3 

k ( t - T ) a 2 = E - e sin E: , 


( 2 ) 


where a is the mean radius, e the eccentricity, and T the time at which the eccentric 
anomaly vanishes. The relation between eccentric anomaly and true anomaly is given by 


e - cos E 

cos v — — , 

1 - e cos E 


where v is the true anomaly. 

Gauss devised a method for solving these transcendental equations for the orbit 
parameters in terms of two points on the ellipse and the time difference between them. 
The method developed by Gauss was extended to other oibit classes by Olbers. Let r , , 
&i and r 2 ,p 2 be two positions of the satellite at times i, and t 2 , respectively; r is the 
distance from the center of the earth, and p is the angular position of the satellite in the 
plane of the orbit measured with respect to the line of ncdes. In terms of these quantities 
the Gauss solution gives the following expression for the mean radius: 

rj + r 2 - 2 • r j r 2 cos g co 5 f 

a = 

2 sin 2 g 

The eccentricity is given by 


if 


n 
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In this equation p is the semilatus rectum 


nr. r, sin 2f \ 2 

P = \ ^2 - ‘i) J 

where f = (fi 2 -/3 1 )/2, k - >mg, G is the gravitational constant, and M is the mass of the earth. 
The quantities r and - are calculated as follows: 


sin 


R 
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where 


l = 


'l + r 2 


4 v y r 1 r 2 cos f 


1 

2 ’ 


k( t 2 t < ) 


^2 v r j r 2 cos f j 




and r; is the root of the following equation: 


T r “ 6 7 i - hr; = 0 . 

g 


The quantity h in this equation is given by 


6 + 1 + 


where 


£ 
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with 



An iterative method was used to solve these equations. 


APPROXIMATION FOR SMALL ECCENTRICITY 

If the eccentricity is less than 0.01, it is sometimes convenient to replace the 
preceding equations by new ones based on the assumption that a circular orbit is a good 
approximation for the motion. The development can be carried out either by expanding 
the Gauss-Olbers solution in powers of terms of the order of the eccentricity, or by 
expanding the original Equations 1 and 2 in the same manner. The results of this devel- 
opment are as follows: 
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The mean motion is calculated from 


N = N 0 



r 2 + r l 


where 


N 0 


2f 


t 


1 


and 

r “i - 

k(t 2 - t t ) 3 

a 0 " [_ 2f 

The eccentricity is then obtained from 


in which 


e 



CSC 



+ tan 2 



2(a - fl ) 
r 2 - r l 


1 
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g 2 - Bl ' N(t 2 “ l l)- 


GEOMETRICAL METHOD 

In the geometrical method the working equations are much simpler, although they 
yield less accurate results. The purely geometrical slatement that three points and the 
focus completely determine an ellipse is used here. H x and y are the Cartesian 
coordinates of the points of the ellipse in the plane of the orbit, then these must satisfy 
the following equation: 

/(x + 2ae cos £) 2 + (y + 2ae sin £) 2 : 2a - >/x 2 + y 2 , 

which constitutes a statement of the basic property of he ellipse. Here 


with 


/x 2 + y 2 + bx + dy - R : 0 , 
b = e cos 8 , d = e sin 8 , R = a( 1 - e 2 ) , 


(3) 


where e is the eccentricity, 8 the argument of perigee, and a the mean radius. 

Equation 3 is linear in the three unknowns b, d, and R, and may be solved directly 
if three points (x, y ) are given by the radar observations. If there is a redundancy of 
data points, the unknowns may be fitted by least squares or by direct averages. 
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Appendix A 

REDUCTION OF RADAR DATA 

Let Sj , Ej, denote slant range, elevation, and bearing at time t 4 . Let R 0 be 
the distance from the center of the earth to the radar station of latitude <p 0 and longitude 
\ 0 . Then, the distance of the satellite from the center of the earth is 

r i 2 = Si 2 + R 0 2 + 2R 0 S ; sin E. . 

The latitude *. and longitude x. of the satellite are computed by means of 

sin 0 £ = sin 0 O cos Q i + cos <p Q sin Qj cos , 

sin (\. - X Q ) sin B i 
sin Q i cos <p. 

Next, transfer to a set of ordinary spherical coordinates, correcting for the earth's 
rotation: 

K [ = K i + " t O ’ 

• 77 

" 2 " * 

Here is the angular velocity of the earth, and t k is the time of the particular 
observation at which the line of nodes is computed. Then the Cartesian coordinates in the 
fixed frame are given by 

X L = r . sin 0. cos , 

Y. -- r | sin c p! sin , 

t 

Z. = r. cos 0. „ 

i i 1 

where the x-axis passes through the Greenwich meridian at time t k and the z-axis is 
the polar axis. The plane of the orbit is found by a least-squares fitting. 

Next, coordinates in the plane of the orbit are calculated by means of 

Xj = rj cos /3. , 

Y = i-j sin /3. , 

COS 0. 

where sin /3. = a ‘ - and a is the inclination angle. 

The planar elements of the ellipse are now computed by using the Cartesian coordi- 
nates Xj, Yj for the geometrical method and the polar coordinates rj, p i for the dynamical 

method. 
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